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1. PROJECT DESCRIPTION 

The funded project seeked a unified de- 
scription of the ionization, physical structure, 
and evolution of the inter galactic medium 
(IGM) and quasar intervening absorption 
systems. We proposed to conduct theoretical 
studies of the IGM and QSO absorbers in the 
context of current theories of galaxy forma- 
tion, developing and using numerical and an- 
alytical techniques aimed at a detailed mod- 
eling of cosmological radiative transfer, gas 
dynamics, and thermal and ionization evolu- 
tion. The ionization history of the IGM has 
important implications for the metagalactic 
UV background, intergalactic helium absorp- 
tion, 21-cm tomography, metal absorption 
systems, fluctuations in the microwave back- 
ground, and the cosmic rate of structure and 
star formation. 

All the original objectives of our program 
have been achieved, and the results widely 
used and quoted by the community. Indeed, 
they remain relevant as the level and com- 
plexity of research in this area has increased 
substantially since our proposal was submit- 
ted, due to new discoveries on galaxy forma- 
tion and evolution, a flood of high-quality 
data on the distant universe, new theoret- 
ical ideas and direct numerical simulations 
of structure formation in hierarchical cluster- 
ing theories. The dawning of the epoch of 
galaxies, the physics of the interplay between 
early galaxy formation and the intergalactic 
medium, and possible observational probes of 
the end of the “dark ages” remain some of the 
most debated topics in observational cosmol- 
ogy. The duration of the program was ex- 
tended beyond the original three years as the 
PI took a one-year sabbatical in the UK and 
afterwards moved to the University of Cali- 
fornia Santa Cruz. A graduate student, C. 
Porciani, has been funded by this proposal. 
Below is a list of some of the most recent re- 


sults obtained. 

2 PHYSICAL PROPERTIES OF THE LYMAN 
ALPHA FOREST 

Avery Meiksin and Mike Norman studied 
the origin and physical nature of the Lya for- 
est absorption systems as found in hydrody- 
namics! simulations of the IGM in a stan- 
dard Cold Dark Matter cosmology (Zhang 
et a!. 1998). The structures of the systems 
that give rise to the Lya forest span a wide 
range in morphologies. At the cosmological 
average density, the characteristic morphol- 
ogy is cell-like with underdense regions sepa- 
rated by overdense sheet-like partitions. The 
lowest density contours tend to enclose amor- 
phous, isolated regions. They found that the 
principal structures of the IGM are in place 
by z = 5, with the evolution in the IGM ab- 
sorption properties due primarily to the ex- 
pansion of the universe and the changing in- 
tensity of the photoionizing background radi- 
ation field. Virtually all the baryons in the 
simulation fragment into structures that we 
can identify with discrete absorption lines, 
with at most 5% remaining in a smoothly dis- 
tributed component (the Gunn-Peterson ef- 
fect). Given reasonable assumption on the 
spectrum of the ionizing UV background, the 
simulations were able to reproduce the statis- 
tics of the carbon and silicon measurements 
of the Lya forest at z 3 and the intergalac- 
tic He II opacity measurements of HST and 
HUT 

3. THE TEMPERATURE OF THE IGM 

Numerical N-body /hydrodynamics simula- 
tions of structure formation in the IGM, 
within the framework of CDM-dominated 
cosmologies have provided a definite picture 
for the origin of the Lya forest, one of an 
interconnected network of sheets and fila- 
ments with virialized systems (halos) located 
at their points of intersection. Most simula- 
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tions to date assume that the IGM is pho- 
toionized and photoheated by a UV radia- 
tion background close to that inferred from 
quasars (Madau, Haardt, & Rees 1999). The 
width of the lines, as measured by the b pa- 
rameter of a Voigt profile, is set by thermal 
broadening, peculiar velocities, and Hubble 
expansion across the filaments. While the 
first simulations showed good agreement with 
the observed line statistics, more detailed 
studies at higher numerical resolution have 
revealed a serious conflict with the data: the 
models predict median values for the b dis- 
tribution of about 20 kms -1 at z = 3, com- 
pared with the observed median of 30 km s _1 . 
One possible way to restore the agreement 
with the observations is to increase the tem- 
perature of the IGM, adding to the thermal 
broadening. We pointed out (Madau & Ef- 
stathiou 1999) that Compton heating of in- 
tergalactic gas by the hard X-ray background 
(XRB) may help to resolve this discrepancy 
in underdense regions. The rate of gain in 
thermal energy by Compton scattering will 
dominate over the energy input from hydro- 
gen photoionization if the XRB energy den- 
sity is ~ 0.2x/(e) times higher than the en- 
ergy density of the UV background at a given 
epoch, where x is the hydrogen neutral frac- 
tion in units of 10 -6 and (e) is the mean 
X-ray photon energy in units of m e c 2 . The 
numerical integration of the time-dependent 
rate equations showed, however, that in pop- 
ular models for the redshift evolution of the 
extragalactic background radiation this effect 
was not large enough to match the observa- 
tions. Perhaps a more promising mechanism 
that could increase the line-width involves ra- 
diative transfer effects during the (late) reion- 
ization of helium. As discussed below , a ma- 
jor effort has been devoted to develop a radia- 
tive transfer code for treating thermal effects 
during the inhomogeneous reionization. 

4. SIMULATING INHOMOGENEOUS 
REIONIZATION 

A key requirement for numerical simula- 
tions of inhomogeneous reionization is an ef- 
ficient algorithm for computing the radiation 
field J v { x) in three spatial dimensions given 
an arbitrary distribution of absorbers and 
sources, both point and diffuse (e.g., recombi- 


nation radiation). In principle, this is accom- 
plished by solving the equation of radiative 
transfer for the monochromatic intensity 
However, since /„ is a function of seven vari- 
ables (t, x, fi, v), this is impractical on current 
computers if the intensity needs to be evalu- 
ated many times over the course of a hydro- 
dynamical calculation. For this reason, we 
have been exploring physically motivated ap- 
proximations to the transfer equation which 
retain the essential physics while being nu- 
merically economical. 

Norman, Paschos & Abel (1998) developed 
a hydrid method that decomposes the radia- 
tion field into point source and diffuse com- 
ponents. The point source radiation field is 
attenuated along radial rays, while the dif- 
fuse radiation field is computed by solving the 
zeroth and first radiation moment equations 
field with variable Eddington factor closure. 
In the limit of small (compared to the hori- 
zon) volumes, the two moment equations can 
be simplified and combined into a single non- 
linear elliptic equation for the diffuse field. 
The elliptic equation can be solved using a 
multigrid algorithm. 

In Abel, Norman, & Madau (1999) a nu- 
merical method that solves the radiative 
transfer around point sources within a three 
dimensional cartesian grid was developed. 
The method is energy conserving indepen- 
dently of resolution: this ensures the correct 
propagation speeds of ionization fronts. We 
described the details of the algorithm, and 
computed as first numerical application the 
ionized region surrounding a miniquasar in a 
cosmological density field at z = 7. 

5 RADIO SIGNATURES OF HI AT HIGH 
REDSHIFT 

The emission of 21-cm radiation from a 
neutral intergalactic medium (IGM) at high 
redshifts has been discussed in a paper (Tozzi 
et al. 2000) in connection with the thermal 
and ionization history of the universe. The 
physical mechanisms that make such radia- 
tion detectable against the cosmic microwave 
background include Lyman-alpha coupling of 
the hydrogen spin temperature to the kinetic 
temperature of the gas and preheating of 
the IGM by the first generation of stars and 
quasars. Three different signatures were in- 
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vestigated in detail: (a) the fluctuations in 
the redshifted 21-cm emission induced by the 
gas density inhomogeneities that develop at 
early times in cold dark matter (CDM) dom- 
inated cosmologies; (b) the sharp absorption 
feature in the radio sky due to the rapid rise 
of the Lyman-alpha continuum background 
that marks the birth of the first UV sources 
in the universe; and (c) the 21-cm emission 
and absorption shells that are generated on 
several Mpc scales around the first bright 
quasars. Future radio observations with the 
Giant Metrewave Radio Telescope and the 
Square Kilometer Array may shed light on 
the power spectrum of density fluctuations 
at z > 7, and map the end of the “dark age”, 
i.e. the transition from a neutral universe to a 
ionized one populated with radiation sources. 

6. PHOTON CONSUMPTION IN MINIHALOS 
DURING REIONIZATION 

The process of reionization and its im- 
pact on several key cosmological issues (from 
the role it plays in allowing collapsed ob- 
jects to cool and make stars, to determin- 
ing the small-scale structure in the temper- 
ature fluctuations of the CMB) has recently 
received much theoretical attention, and has 
been studied by several authors using semi- 
analytic models and three-dimensional nu- 
merical simulations. In general, these works 
aim to follow the time evolution of the fill- 
ing factor of ionized (H ii ) regions, based on 
some input prescriptions for the emissivity 
and spectra of the ionizing sources. 

More recent works have focused on the in- 
creased rate of recombinations in a clumpy 
medium relative to a homogeneous one. 
These studies have left significant uncertain- 
ties on the details of how reionization pro- 
ceeds in an inhomogeneous medium. Since 
the ionizing sources are embedded in dense 
regions, one might expect that these dense 
regions are ionized first, before the radiation 
escapes to ionize the low-density IGM. Alter- 
natively, most of the radiation might escape 
from the local, dense regions along low col- 
umn density lines of sight. In this case, the 
underdense ‘voids’ are ionized first, with the 
ionization of the denser filaments and halos 
lagging behind. 

At the earliest epochs of structure for- 


mation in CDM cosmologies the smallest 
nonlinear objects are the numerous small 
halos that condense with virial tempera- 
tures between the cosmological Jean temper- 
ature and ~ 10 4 K. Such “minihalos” are 
not yet fully resolved nor is the process 
of their photoevaporation captured in large- 
scale three-dimensional cosmological simula- 
tions. In Haiman et al. (2001) we have used 
semi-analytic methods combined with three- 
dimensional numerical simulations of individ- 
ual photoevaporating minihalos to (1) quan- 
tify the importance of high-redshift miniha- 
los as sinks of ionizing photons; and (2) as- 
sess whether by extrapolating to early times 
the known population of galaxies and quasars 
a sufficient number of UV photons are pro- 
duced for hydrogen reionization. If reion- 
ization occurs at sufficiently high redshifts 
( z r > 20), the intergalactic medium is 

heated to ~ 10 4 K and most minihalos never 
form. On the other hand, if z T 20 as it 
appears more likely, then a significant frac- 
tion (,> 10%) of all baryons have already col- 
lapsed into minihalos, and are subsequently 
removed from the halos by photoevaporation 
as the ionizing background flux builds up. 
This process requires a significant budget of 
ionizing photons, about 10-20 of them per 
background hydrogen atom. This exceeds the 
production by a straightforward extrapola- 
tion back in time of known quasar and galaxy 
populations by a factor of up to ~ 10 and ~ 3, 
respectively (Haiman et al. 2001). 

7. THE EARLIEST LUMINOUS SOURCES 

Prior to complete reionization at redshift 
z r , sources of ultraviolet radiation will be 
seen behind a large column of intervening 
gas that is still neutral. In this case, be- 
cause of scattering off the line-of-sight due 
to the diffuse neutral IGM, the spectrum of 
a source at z em > z r should show the red 
damping wing of the Gunn-Peterson absorp- 
tion trough at wavelengths longer than the 
local Lya resonance, A c bs > \»(1 + z em ) i 
where X a — c/v Q — 1216 A. At z e m 6, this 
characteristic feature extends for more than 
1500 kins -1 to the red of the resonance, and 
may significantly suppress the Lya emission 
line in the spectra of the first generation of 
objects in the universe. Measuring the shape 
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of the absorption profile of the damping wing 
could provide a determination of the density 
of the neutral IGM near the source. 

In Madau & Rees (2000) we focused on the 
width of the red damping wing - related to 
the expected strength of the Lya emission 
line - in the spectra of very distant QSOs 
as a flag of the observation of the IGM be- 
fore reionization. We assessed, in particu- 
lar, the impact of the photoionized, Mpc-size 
regions which will surround individual lumi- 
nous sources of Lyman continuum radiation 
on the transmission of photons redward of the 
Lya resonance, and shown that the damp- 
ing wing of the Gunn-Peterson trough may 
nearly completely disappear because of the 
lack of neutral hydrogen in the vicinity of 
a bright QSO. The effect of this local pho- 
toionization is to greatly reduce the scatter- 
ing opacity between the redshift of the quasar 
and the boundary of its H n region. The de- 
tection of a strong Lya emission line in the 
spectra of bright QSOs shining for ,> 10 7 yr 
cannot then be used, by itself, as a constraint 
on the reionization epoch. The first signs 
of an object radiating prior to the transition 
from a neutral to an ionized universe may be 
best searched for in the spectra of luminous 
sources with a small escape fraction of Lyman 
continuum photons into the IGM, or sources 
with a short duty cycle (like, e.g., gamma-ray 
bursts) . 

8. EARLY METAL ENRICHMENT OF THE IGM 

Understanding the origin of the chemical 
elements, following the increase in their abun- 
dances with cosmic time, and uncovering 
the processes responsible for distributing the 
products of stellar nucleosynthesis over very 
large distances are all key aspects of the evo- 
lution of gaseous matter in the universe. One 
of the major discoveries with Keck concern- 
ing the IGM has been the identification of 
metal absorption lines associated with many 
of the Lya forest systems at z = 3 — 3.5. 

Prom a theoretical perspective, it is un- 
known whether the existence of heavy ele- 
ments in the low-density IGM points to an 
early (z > 6) enrichment epoch by low-mass 
subgalactic systems (Madau, Ferrara, & Rees 
2001), or is rather due to late pollution by the 
population of star-forming galaxies galaxies 


known to be already in place at z = 3. Mas- 
sive stars will deposit both radiative and me- 
chanical energy into the interstellar medium 
of protogalaxies. A complex network of ‘feed- 
back’ mechanisms is likely at work in these 
systems, as the gas in shallow potential is 
more easily blown away, thereby quenching 
star formation. To start addressing some 
of these issues and simulate the process of 
blow-away, we have developed (Mori, Fer- 
rara, & Madau 2002) a 3D Eulerian code 
that solves the hydrodynamic equations for 
a perfect fluid in Cartesian geometry. To 
deal with very different length scales in our 
simulation we have adopted a ‘nested grid 
method’ with six levels of fixed Cartesian 
grids. The scheme has a wide dynamic range 
in the space dimension. The results of our 
simulations show that SN driven pregalac- 
tic outflows may be an efficient mechanism 
for spreading metals around. Furthermore, 
as the blastwaves produced by supernova ex- 
plosions reheat the surrounding intergalactic 
gas and enrich it with newly formed heavy 
elements, they appear to inhibit the forma- 
tion of surrounding low-mass galaxies due to 
‘baryonic stripping’ (Scannapieco, Ferrara, &: 
Madau 2002). 
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